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ABSTRACT 

A model is developed for the time dependent electromagnetic - radio to gamma-ray - emission of active 
galactic nuclei, specifically, the blazars, based on the acceleration and creation of leptons at a propa- 
gating discontinuity or front of a self-collimated Poynting flux jet. The front corresponds to a discrete 
relativistic jet component as observed with very-long-baseline- interferometry (VLBI). Equations are 
derived for the number, momentum, and energy of particles in the front taking into account syn- 
chrotron, synchrotron-self-Compton (SSC), and inverse-Compton processes as well as photon-photon 
pair production. The apparent synchrotron, SSC, and inverse Compton luminosities as functions of 
time are determined. Predictions of the model are compared with observations in the gamma, optical, 
and radio bands. The delay between the high-energy gamma-ray flare and the onset of the radio is 
explained by self- absorption and/or free- free absorption by external plasma. Two types of gamma-ray 
flares are predicted, Compton dominated or SSC dominated, depending on the initial parameters in 
the front. The theory is applied to the recently observed gamma-ray flare of the blazar PKS 1622-297 
(Mattox et al. 1996). Approximate agreement of theoretical and observed light curves is obtained for 
a viewing angle 9 i, s ~ 0.1 rad, a black hole mass M ~ 3 x 1O 9 M0, and a magnetic field at the base 
of the jet B ~ 10 3 G. 

Subject headings: active — galaxies: quasars — galaxies: jets — galaxies: gamma-rays — galaxies 



1. INTRODUCTION 

New understanding of the nature of Active Galactic 
Nuclei (AGNs) has come from the discovery of the high 
energy gamma-ray radiation in the range 50 — 10 4 McV 
by the EGRET instrument on the Compton gamma-ray 
Observatory (e.g., Hartman et al. 1992; Mattox et al. 
1993; Thompson et al. 1993). This radiation is observed 
from a sub-class of AGNs termed blazars, which include 
Optically Violently Variable (OVV) quasars and BL Lac 
objects and which show strong variability in all wave- 
bands from radio to gamma. Many of the objects reveal 
'super-luminal' jets in VLBI maps which indicate that 
we observe matter of the jet pointed nearly towards us 
and that the jet matter moves with relativistic speed. 

Prior to the Compton Observatory measurements, pre- 
diction of strong, collimated gamma-ray emission from 
AGN relativistic jet sources was made by the electromag- 
netic cascade model of Lovelace, MacAuslan, and Burns 
(1979) and Burns and Lovelace (1982). More recently, a 
number of theoretical models have been developed to ex- 
plain the observed gamma-ray emission of AGNs (see re- 
view by Sikora 1994). In most of the models the gamma- 
ray radiation is ascribed to inverse-Compton scattering 
of relativisitic electrons and possibly positrons (Lorentz 
factors 7 ~ 10 2 — 10 3 ) of a jet having relativistic bulk 
motion (Lorentz factor T ~ 10) with soft photons (en- 
ergies ~ 1 — 10 2 eV). The soft photons can arise from 
the synchrotron emission of the relativistic electrons in 
the jet as in the synchrotron-self-Compton (SSC) mod- 
els (Maraschi, Ghisellini, and Celotti 1992; Marscher and 
Bloom, 1992) or from the direct or scattered thermal ra- 



diation from an accretion disk (Dermer, Schlickeiser, and 
Mastichiadis 1992; Blandford 1993; and Sikora, Begel- 
man, and Rees 1994), or from a single cloud (Ghisellini 
and Madau 1996). The spectrum of blazars in a pair 
cascade model was calculated by Levinson and Bland- 
ford (1995) and Levinson (1996). In different models, 
ultra high energy protons (Lorentz factors > 10 6 ) are 
postulated to cause a cascade, the product particles of 
which produce the observed radiation (Mannheim and 
Biermann 1992), or the jets considered to be ultra rela- 
tivistic with bulk Lorentz factors T > 10 4 (Coppi, Kartje, 
and Konigl 1993). 

Here, we propose that the main driving force for the ob- 
served superluminal jet components is a finite amplitude 
discontinuity in a Poynting flux jet. A brief discussion 
of the model was given earlier (Lovelace and Romanova 
1996). A rapid change in the Poynting jet outflow from a 
disk can result from implosive accretion in a disk with an 
ordered magnetic field (Lovelace, Romanova, and New- 
man 1994). Propagation of newly expelled EM field and 
matter from the disk with higher velocity than the old 
jet can under certain conditions lead to the formation 
of a pair of shock waves as in the case of non-relativistic 
hydrodynamic flow (Raga et al. 1990). Particle accelera- 
tion in the front may result from the shocks or from anni- 
hilation and/or reconnection of oppositely directed mag- 
netic fields in the front (Romanova and Lovelace 1992 
and Lovelace, Newman, and Romanova 1996). 

Section 2 derives a complete set of equations for the 
front, specifically, equations for the total number of par- 
ticles, the total momentum, energy, and magnetic flux 
in the front. Section 3 discusses and summarizes the be- 



2 



havior of the solutions of the front equations. Section 4 
gives comparisons of predictions of the theory with ob- 
servations. Conclusions of this work are given in Section 
5. 

2. THEORY 

We use an inertial, cylindrical coordinate system 
(r, <f>, z) with the origin at the black hole's center and 
the axis normal to the accretion disk as shown in Fig- 
ure 1. This is referred to as the 'lab frame'. We consider 
that Poynting flux jets propagate symmetrically away 
from the disk in the ±z directions, but focus our at- 
tention on the approaching +z jet. A Poynting flux jet 
is self-collimated, with energy, momentum, and angular 
momentum transported mainly by the electromagnetic 
fields (Lovelace, Wang, and Sulkanen 1987). The colli- 
mation results from the toroidal magnetic field at the 
edge of the jet. A steady Poynting jet can be char- 
acterized in the lab frame by its asymptotic (z 3> r Q ) 
magnetic field = — B[r a /rj(z)], and electric field 
E r = —(v/c)B{r /rj(z)] at the jet's edge, r — rj(z), 
where r a is the jet's radius at z = 0, v = const, is the 
jet's axial velocity, and B is the lab frame field strength 
at z = 0. 

The jet plasma consists of both ions and leptons (elec- 
trons and positrons) with the ratio of leptons to ions fu 
The initial jet radius is taken to be r a = 6GM/c 2 , where 
M is the black hole mass. The energy flux (luminos- 
ity) of the +z jet is the Poynting flux Lj = vB 2 r 2 /8 = 
constant. 

We propose that the Poynting jet from the disk changes 
abruptly at time t = 0. That is, the jet parameters 
change from values with subscript (1) to subscript (2) 
at t = 0, ( Bi, na, vi, fin) — » ( B 2 , n i2 , v 2 , fm) ■ In 
actuality the change will be with a time scale determined 
by the disk dynamics as in the implosive accretion model 
of Lovelace et al. (1994). This time scale may be as short 
as r Q /c in the lab frame. In the present work we consider 
b = Bi/B 2 < 1, v = nn/n i2 , Vi < v 2 , f H1 = f U2 , 
where m is the number density of ions. The change in 
the jet parameters produces a 'front' which propagates 
outward as indicated in Figure 1 . The front may involve 
a pair of shock waves, one for the incoming old jet matter 
and the other for the incoming new jet matter. We let 
Z(t) = z(t)/r denote the dimensionless axial distance 
of the front. We also use lab time T — tj (r /c), speed of 
the front U(T) = dZ/dT = (dz/dt)/c, and bulk Lorentz 
factor r = (1 - U 2 )-^. We let Vi, 2 = vi, 2 /c and ri, 2 = 
(1 — V 2 2 )~z. The time measured in the front frame is 
V = t'/(r /c) with dT' = dT/Y{T). The initial values 
at T = V = are Z = and U = V x . We also use 
the time T" = t" I (r /c) measured by a distant observer 
oriented at an angle 8 b s < ir/2 to the z— axis; dT" — 
dT[l - U cos{9 obs )}. Further, the observed time T'" for 
a cosmological source is given by T'" = T"(l + z) with z 
the redshift. 

2.1. Number of Particles 
The continuity equation for ions in the front frame is 
dn'j _ d(nW) 

dt' dz' ' 

where n\, v', and t' are the number density, velocity, and 
time in the front frame. We integrate this equation over a 
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cylindrical 'pill box' of radius r > rj(z) and axial length 
Az' and use the Lorentz transformations to obtain an 
equation for the total number of ions in the front (T) , 



dT 



-^ [AF 2 fl r (Ay 2 ) + ^Ay 1 if(AI/ 1 )]. (2) 



Here, Ni = irr 2 rii 2 (z)r — nrg(m2) z=o is a constant 
by our earlier assumptions. Also, H(x) = 1 for x > 
and H(x) — for x < 0. If Az' were a constant, then we 
would have AVi = U—Vi and AV 2 = V 2 — U, where upper 
case quantities are dimensionless. However, the plasma 
in the front expands freely with speed C' s (normalized by 

c) so that AZ' = AZ' +2 J^' dT"C' s {T"), and this results 
in modified expressions for AVi j2 . In the radial direction 
the plasma also expands freely with sound speed C' s . C' s 
is the compressional wave speed in the front frame. For 
the cases studied, C' s is determined essentially by the 
particles in that the particle energy density is larger than 
that of the magnetic field. If both the ions and leptons 
are highly relativistic in the front frame, C' s = 1/V3, 
whereas for non-relativistic ions but relativistic leptons 
C' s < 1/ \/3- Magnetic pinching of the front is estimated 
to be small. Note that the number density of ions in the 
front frame is simply n\ = Ni/ (nr^R 2 AZ'). The electric 
field in the front frame |E'| is small compared with |B'|. 

If fin = fli2 and there is no e± pair production in the 
front, then the total number of leptons in the front is sim- 
ply = fuNf. However, in the general case considered 
here, 
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N io [f U2 AV 2 H 2 + fuivAViHi] + 1 ^ 

where the H's are the same as in equation (2). The 
main contribution to the pair production for the con- 
ditions considered is from collisions of synchrotron and 
SSC photons, 

\ST> 
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de 2 
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For rough estimates, we can write ei = ^ 2 hw' , where 



uj' = e\B'\/(m e c) is the cyclotron frequency in the front 
frame, and e 2 = 7 2 ei. Thus, a rough condition for 
an electron to give pair production is 7 > ^ pa ir = 

(mc 2 /hLu' )^ w 3.5x 10 3 (10 3 G/|B'|)3. Electron-positron 
recombination is negligible for the conditions considered. 



2.2. Momentum Conservation 
In the front frame, 

-^ = -^ +9 rav + rad, 



(4) 



where Tq z and T' zz are components of the energy momen- 
tum flux density tensor in the front frame, and 'grav' and 
'rad' denote gravitational and radiative force contribu- 
tions not included in T[^ . We integrate (4) over the same 
pill box to obtain 



(Nimai + Nim e ji) r 
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Fig. 1. — Sketch of the geometry of a propagating front in a Poynting flux jet. 



if 



(grav + rad) . 



(5) 



7i and 7; denote averages over the ion and lepton dis- 
tribution functions in the front frame. In this frame the 
distributions are assumed isotropic and that for electrons 
is assumed the same as for positrons. The initial values 
of % and 7; are considered to be close to unity. A Lorentz 
transform gives (T' zz ) g = T 2 (T ZZ - 2UT oz + U 2 T ocl ) s , 
where s — 1,2. The lab frame components of the energy- 
momentum tensor for a Poynting flux jet are 



T — 



8tt 

E r B^ 
4ir 



8tt 



- 2(V) 6 



B 2 ^ 

8W, 



(6) 



where B = at the edge of the jet. Thus, in equation 
(5) we have 

- [T»] = ( [(1 + U 2 ){1 + V 2 )- 4UV 2 ] - 



nr 2 r, 



^[(l + U^il + V^-WVi] 



(7) 



where b 2 = (B1/.B2) 2 , and r 2 B\ = const = rl(B 2 ) 2 z=0 . 
We assume b 2 < 1. The terms in the square brackets in- 
volving V2 represent the push from the new Poynting jet, 
while those involving V\ are for the push in the opposite 
direction from the old Poynting jet. A small modifica- 
tion of equation (7) similar to that of (2) is required to 
account for the free expansion of the plasma of the front. 
Note that the electromagnetic field contribution to the 
momentum of the front is small because |E'| <C |B'|. 

Dividing equation (5) by AM Q = A^ojm; + iV;(0)m e 
gives 



Nimai + Nim e ji \ 3 dU 



AM„ 



r^ = ,r 2 {[(..)] 



AM, 



^2 / (drav + rad) , 



(8) 



where /i = r^(B2)1 =0 / (SAM c 2 ) 3> 1 is a dimensionless 
measure of the strength of the Poynting jet. The brackets 
{} denote the same quantity as in equation (7). The 
driving term cx \i vanishes if both Ti — > T 2 and b — > 1. 
If r, Ti, and T2 are all much larger than unity, it is clear 
from equations (7) and (8) that a steady state is possible 
if b 2 > (ri/r 2 ) 4 (recall that b 2 < 1); that is, the push of 



4 



the old and new jets balance and the Lorentz factor of 
the front is T = Ti[(l - b 2 )/(b 2 - (ryr 2 ) 4 )]3. If b 2 is 
smaller than (Ti/I^) 4 , then no balance is possible and 
the Lorentz factor of the front T increases without limit. 

The gravitational force in equation (8) can be written 
as J grav = -CM (iVjraj + Nim e ) / \r 2 + z 2 ) . The radia- 
tive force depends in general on the geometry and energy 
distribution of the background radiation field of the cen- 
tral region of the AGN. Dermer et al. (1992) consider 
the case where the radiation comes from the disk, while 
Sikora et al. (1994) argue that the radiation field is from 
disk radiation scattered by a distribution of clouds orbit- 
ing the central object. We adopt a rough parameteriza- 
tion of the radiation fields of Dermer et al. and Sikora et 
al. The average photon energy is denoted e p h, the total 
luminosity L p h, and the characteristic radius of the spa- 
tial distribution r p ^ . The force due to this radiation field 
is J rod = (Nio T )L ph V pz T T /[^{r 2 ph + (z - ar ph ) 2 )c], 
where a = for the Dermer et al. model and a = 1 
for that of Sikora et al. Here, T> pz = Y 2 [cos(9 p h) — 
U] \cos(9 p h) — U\ is the Doppler factor which accounts 
for the change in the flux and the change in the z mo- 
mentum of the photons between the lab and front frames, 
cos(9 ph ) = (z - ar ph )(r 2 h + (z - ar ph ) 2 )-2 . u T denotes 
the Thomson cross section. The contribution of leptons 
with 7 > m e c 2 /e p h to (Ni<jt) in the radiative force is 
reduced in that the Klein-Nishina cross section applies. 
Tt = (1 — e~ TT )/TT with tt = n\ r a R o~t the Thomson 
optical depth of the front. For the conditions considered 
here, tt < 1. For small axial distances, cos(9 p h) < U, 
the radiative force acts as a drag whereas at larger dis- 
tances it gives a push. 

2.3. Energy Conservation 
In the front frame, 



dt' 



dV 
2_Z 21 

' dz' 



syn — ssc — Com 



(9) 



where the last three terms represent the energy loss rates 
due to synchrotron radiation, inverse-Compton scatter- 
ing off of synchrotron photons, and inverse-Compton 
scattering off of the above mentioned background pho- 
tons. Following our previous method, we integrate (9) 
over the volume of the front to obtain 



dT 



N^l-^mi + N^l 



l)m e + ^ 
c z 



7rr 2 r, 



~ [ T 0z] ~~^jf ( Syn + SSC + Com) 



(10) 



where W' B = (rl/Ac 2 )R 2 AZ'(B') 2 is the magnetic field 
energy in the front. The magnetic field B' is discussed 
below in subsection 2.4. The driving term oc u vanishes 
if both ri — > T 2 and b — > 1. Because the acceleration 
process (es) in the front is not known, we make the well- 
known supposition (for example, Pacholczyk 1970) that 
the kinetic energy in the ions is a constant factor k times 
that in the leptons, Ni(ji — l)mt — kNi(ji — l)m e . As 
a result, ji is determined in terms of 7/. Therefore, in 
the following 7 without a subscript refers always to the 
lepton Lorentz factor. 



A Lorentz transform gives T' oz = T 2 [(l + U 2 )T oz - 
UT QO - UT ZZ ] , so that 

2 R 2 r 2 e 

J [(l + U 2 )V 2 -U{l + V 2 )]H 2 



-Kr 2 r n 



c 3 



rpl 



4c 2 



-6 2 [(l + i7 2 )F 1 -t/(l + F 1 2 )]tf 1 J . (11) 



Dividing equation (10) by AM a gives 



d 



AM dT 



l)(l + fc)m e + — s 



2^r 2 { [(..)]} 



(syn+ ..) , 



(12) 



AM„c 3 

where the brackets {} denote the same quantity as in 
equation (11). The driving term oc u vanishes if both 
Ti -> T 2 and 6 -» 1. 
The synchrotron energy loss rate in the front frame is: 



/ 



/\2 



syn = -g-^JV,( 7 ; ) 2 c(^ 



0* 



syni 



(13a) 



where r P 



e 2 /(m e c 2 ) is the classical electron radius, 
and ( 7 ) 2 = J™ c?77 2 /;/ Ni, where 71 is the Lorentz 
factor below which synchrotron self-absorption becomes 



strong as discussed below. T s . 



[1 CXp( 1~syn)\l Tsyn 



is synchrotron opacity factor with T syn the optical 
depth for synchrotron photons (Pacholczyk 1970). The 
synchrotron-self-Compton energy loss rate is: 



/ 



4 (AT — f MN^O-T ^iB') 2 ^ 

ssc = -(^a T )(7) 2 c( H— )^ 



(136) 

The energy loss rate in the front frame due to Compton 
scattering off of the background photons is: 

J ^m = r^ar^ciU'^TT , (13c) 

where U' ph = L uv V 2 h l\Kr 2 c(R 2 Yl + (Z - aR ph ) 2 )] is the 
background photon energy density in the front frame, 
and T> p h = T[l — Ucos(9 p h)\ is the Doppler factor for 
photons in the front frame relative to the disk. 

2.4. Magnetic Flux Conservation 

The magnetic field in front frame B' is determined by 
taking into account the influx of magnetic flux with the 
new jet matter (subscript 2) and the old jet matter (sub- 
script 1). We let = RAZ'B' denote the toroidal flux 
(in Gauss) in the front. Then 

^ = B [AV 2 H(AV 2 ) + 6AVi JT(AVi)] , (14) 

where B a = (B 2 ) z= q and the H's are the same as in 
equation (2). If b < 0, there is annihilation of magnetic 
flux in the front (see Romanova and Lovelace 1992; and 
Lovelace, Newman, and Romanova 1996). 

2.5. Lepton Distribution 

The lepton distribution function suggested by obser- 
vations has a hard power law form in the main energy 
containing range, say, /; = ifi/7 2 for 71 < 7 < 72 with 
1 <SC 71 <C 72- For 7 > 72, the distribution is steeper, 
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say, fi = K 2 /j 3 , and at even larger energies, 7 > 73, /; 
is even steeper. For 7 < 71, /; is assumed negligible (see 
Figure 2). Thus, fi is characterized mainly by 71 and 

72. We have 7 = 71 [1 + £72(72/71)], and 7 2 = C7172, with 
C = 1 + ^(73/72). The relevant value of 71 is the Lorentz 
factor below which synchrotron self-absorption becomes 
strong. Self-absorption as treated for example by Pa- 
cholczyk (1970) then gives 71 = ji(B', nj, R). With 71 
known, we have 72 = Jiexp(j/^i — 1). The value of C 
depends only weakly on 73/72 and we take 73/72 =7/71- 

3. RESULTS 

Equations (2), (3), (8), (12) and (14) have been 
solved numerically to obtain the time dependences of 
the physical variables. Here, we first discuss the time- 
dependences for a reference case where the black hole 
mass M = 10 9 M Q so that r Q = 8.9 x 10 14 cm, /x = 
r 3 o (B 2 ) 2 z=0 /(8AM o c 2 ) = 15, B = 10 3 G, b = \B l /B 2 \ = 
0.5, Ti = 8, T 2 = 18, v = n a /n i2 = 0.44, k = 1, and 
the initial lepton-ion ratio / z ° = 5. For the background 
photons, L uv = 10 46 erg/s, e ph = 10 eV, r ph — 0.1 pc, 
and a = 0. The angle between the observer and the jet 
axis is Qobs = 0.1 rad. All results arc presented for the 
dependence on time (i") seen by a distant observer where 
it is recalled that dt" = dt[l — Ucos(9 bs)]- The redshift 
dependences are not included in the equations. 

3.1. Time Dependences 

The velocity of the front U, initially equal to Vi, in- 
creases rapidly and approaches the equilibrium value Uoo 
which corresponds to the bulk Lorentz factor — 
10. This is shown in Figure 3. The apparent veloc- 
ity of the front (with respect to the nucleus) V app — 
Usin(9 b s )/[1 — Ucos(6 b s )] approaches 9.9. The Dopplcr 
boost factor 5 = 1/[T(1 — Ucos{9 b s ))\ approaches 10. 

The total number of ions in the front grows with time, 
because under most conditions the front moves faster 
than the old jet matter and slower than the new jet mat- 
ter, thus accumulating particles from both sides. This 
is shown in Figure 4. For the case shown there is no 
pair production. Thus, the ratio of leptons to ions / ; ° 
decreases with time and approaches unity. Other cases 
were calculated where pairs are created and the number 
of leptons grows so that fu becomes much larger than its 
intial value. 

Although the total number of particles in the front 
grows with time, the number densities decrease with time 
due to the expansion of the volume of the front. This is 
shown in Figure 5. The number densities of synchrotron, 
SSC, background, and inverse-Compton photons in the 
front also decrease with time as shown in Figure 5. 

Figure 6 shows the evolution of the different character- 
istic Lorentz factors, 71, 72, and 73, of the lepton distribu- 
tion. Initially, the leptons are assumed to have the same 
small energy corresponding to a Lorentz factor 7 = 1.01. 
Later, when the average lepton Lorentz factor (7) be- 
comes larger than that corresponding to self-absorption 
(71), a complete spectrum forms extending from 71 to 

73. Initially, 71,72, and 73 grow rapidly and later more 
slowly. The top curve in Figure 6, j pa i r , corresponds to 
the energy of leptons at which pair creation sets in. 

The total apparent luminosities due to synchrotron, 
SSC, and Compton processes for a distant observer at an 



angle 9 b s to the line of sight are calculated including the 
solid angle boost factor S 2 between the front frame and 
the observer, where S = [T(l — Ucos(9 b s ))]^ 1 ■ These are 
shown in Figure 7. At early times in an outburst, the SSC 
radiation predominates because the magnetic field in the 
front is strong and the number of synchrotron photons is 
higher than that of background photons. Later, the mag- 
netic field in the front decreases and thus the background 
photons predominate. At the same time, the leptons 
are accelerated in the front and consequently inverse- 
Compton radiation grows, and then later decreases. The 
time of the peak of the inverse-Compton emission de- 
pends approximately linearly on the radius r p h of the 
background photon distribution. Synchrotron radiation 
is at much lower luminosity during the first part of the 
front propagation. The dependences of the maximum lu- 
minosities of the different processes are discussed later. 

The observed frequency is boosted by a factor 5. The 
frequencies of radiation corresponding to different pro- 
cesses and to the characteristic values 71, 72, and 73 of 
the lepton distribution, are shown in Figure 8. One can 
see that the radiation of the front covers most of the fre- 
quency range from the self-absorption limit (the lowest 
curve), up to high gamma ray energies. In many bands of 
the spectrum, there are contributions from two processes. 
For example, in the Egret band (from 10 22 5 Hz to 10 24 5 
Hz) we typically have contributions from both inverse- 
Compton and SSC radiation. Initially, synchrotron radi- 
ation contributes to the band from IR to UV, SSC cov- 
ers the range of frequencies from soft X-ray to EGRET 
GeV energies, and the inverse-Compton covers the range 
from the high-energy X-ray to the EGRET GeV energies. 
Thus, during an EGRET high-energy flare one may ex- 
pect almost simultaneous flashes in wavebands from IR 
to very high energies, excluding the radio band which is 
initially self-absorbed. 

3.2. Dependences of Maxima of Flares 

The maxima of the total apparent luminosities due 
to inverse-Compton (Lc m), synchrotron-sclf-Compton 
(Lssc) and synchrotron processes (L syn ) depend on the 
parameters 9 Q b s , B Q , fj,, M, ffc, r p h, and Luv approxi- 
mately as 



L Com ~ 1.7 x 10 
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if the parameters arc of the order of the normalization 
values. 
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Fig. 2. — Lcpton distribution function in the comoving frame of the front. 71, 72, and 73 are derived from the equations of Section 2. 



The dependences on /j, in equation (15) are approxi- 
mately valid for 10 < \i < 50. For \i < 10, the power is 
higher for SSC and synchrotron radiation, but the same 
for Compton radiation. For fi > 50, the luminosities 
go to a constant. The dependences on B a are valid for 
B a < 2000 G, and become flatter for stronger fields. The 
dependence on / ; ° is about the same as in equation (15) 
for 5 < /;° < 15, and becomes steeper for SSC and syn- 
chrotron radiation for / ; ° > 15 and does not depend on 
/ ; ° for ffo < 4. As for the dependences on o b Sl the L's 
depend more strongly on o t, s than in equation (15) at 
larger angles (8 b s > 0.1) and more weakly for smaller an- 
gles. The dependence on M is approximately the same 
at all parameter ranges. 

Notice that the shape of the light curves are quite dif- 
ferent for different parameters. There are several possi- 
bilities: (1) Both SSC and Compton radiation may be 
in the EGRET band. Then, the shape of the EGRET 
light curve is determined essentially by the SSC radia- 
tion which gives a fast rise of luminosity with typical 
time scale of growth t" = 0.03 days = 43 min (see Fig- 
ure 9a). The rise time may of course be slower due to a 
more gradual transistion from the old jet to the new. (2) 
SSC radiation frequencies may be lower than the EGRET 
band so that only Compton radiation is observed. Then, 
the light curve shape will be determined by Compton ra- 
diation with a maximum at t" = 0.13 days for our refer- 
ence values (see Figure 9b). (3) No radiation may fall in 
the EGRET band. For example, changing [i (with other 



parameters fixed) we get a sharp SSC flare for [i > 35, 
a flatter Compton flare for 20 < yu < 35, and no radi- 
ation in the EGRET band for fi < 20. Thus, for this 
set of parameters, the value [i, which is the initial ra- 
tio of magnetic to rest mass energy-density, should be 
relatively high. The shape of the EGRET light curves 
also depends, for example, on the initial ratio of leptons 
to ions. For 1 < f° { < 5, SSC type flares predominate. 
For 5 < ,/)° < 15, Compton type flares predominate, 
and for / ; ° > 15, no radiation falls in the EGRET band. 
Thus, the shape of the flare may be different in differ- 
ent sources or may change in the same source between 
different outbursts, depending on the magnetic/particle 
rest mass energy-density ratio, or on the initial ratio of 
leptons to ions in the front plasma. 

The maxima of the flares in the SSC and synchrotron 
radiation occur very soon after outburst (much less than 
a day). The maximum of the inverse-Compton flare is 
much later and may be observable, 



O pbs 

0.1 



r p h 
Q.lpc 



(16) 

in days. Note that t" Com has essentially no dependence 
on M. 

Later, after the maxima, all luminosities decrease. The 
SSC luminosity decreases as Lssc ~ (*")~ 2 > the syn- 
chrotron luminosity as L syn ~ (t") -1 , and the inverse- 
Compton luminosity as Lc om ~ (t")~ 4 for t" < 2 — 3 




Fig. 3. — Bulk Lorentz factor T, apparent velocity of the front v apP (normalized to c), and Dopplcr boost factor 5 versus time measured 
by a distant observer for the reference parameters given in Section 3. 
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days, but later for t > 15 days it flattens to ~ (t") . 

The rise of the inverse- Compton luminosity is deter- 
mined by the growth of the number and energy of leptons 
in the front. The front is transparent to inverse-Compton 
radiation from the beginning. At the same time, initially 
plasma of the front is opaque to the synchrotron photons, 
and the rapid growth of synchrotron and SSC radiation is 
determined by the fact that plasma becomes transparent 
to synchrotron photons. The decline of inverse-Compton 
radiation is due to the fact that the energy density of 
background photons decreases in the front frame. An im- 
portant factor in the decrease is the Doppler shift of the 
background photons relative to the front for Z > Rph- 
As for synchrotron and SSC radiation, they decrease with 
time because the magnetic field in the front decreases 
with time. 

4. COMPARISON WITH OBSERVATIONS 

4.1. Gamma-Ray Variability: Comparison with the 
Flare m PKS 1622-297 

Here, we compare theoretical and observed light 
curves. The luminosity in the EGRET band (10 22 5 - 
10 24 ' 5 ) Hz is determined in general by both inverse- 
Compton and SSC processes. We integrated the gamma 
radiation in the EGRET band and got two qualitatively 
different types of light curves (Figures 9a and 9b) . The 
difference is determined by the fact that in some cases 
(depending on parameters) SSC radiation is below the 
EGRET frequency limit and the light curve is deter- 
mined by the inverse-Compton light curve (see Figure 
9b). In the opposite case, when the SSC radiation is 
within the EGRET band, the light curve is determined 
by both processes, but the SSC radiation has a strong 
peak at small time scales, and this determines the shape 
of the light curve (Figure 9a). In the first case, it is more 
probable to observe the fast rise and gradual decrease 
of gamma luminosity, compared with the second case, 
where the rise of luminosity is so fast (maximum reached 
at t" = 0.02 days) that it will be difficult to observe 
and only the 'tail' of the flare will be seen (see, however, 
Mattox et al. 1996). 

EGRET observations of gamma-ray blazars show dif- 
ferent types of flares. In some cases the growth of the 
luminosity is relatively slow compared with the decay 
(Kniffen et al. 1993). In other cases (for example, Hart- 
man et al. 1993), only a decaying light curve is observed 
during the EGRET set of observations. In the third type 
of flare, the rise of the luminosity is much faster than the 
decay (Mattox et al. 1996). 

Time resolved observations of the brightest observed 
gamma-ray flare from blazar PKS 1622-29 (Mattox et al. 
1996) show that during the main flare the flux increases 
by an order of magnitude on a 2 hour time scale, whereas 
the decay is on a 7 hour time scale. Thus the rise is faster 
than the decay. We fitted this flare by a 'Compton type' 
flare using the parameters: M = 3 x 10 9 M Q , B Q = 1700 
G, 9 obs = 0.1 rad = 11.5°, fx = 20, /£ = 6, r ph = 0.3 pc, 
e ph = 10 cV, L uv = 2 x 10 46 ergs/s, and a = 1. This is 
shown in Figure 10. For calculation of the theoretical flux 
we used standard formulae for the luminosity distance at 
the redshift of the source z = 0.815, Hubble's constant 
H Q = 75 km/s/Mpc, and a cosmological parameter q a = 
1/2. Also, we took into account the cosmological redshift 



of the frequencies and the dilation of the time scale of 
the outburst. The value a = 1 of Figure 10 corresponds 
to background photons scattered by clouds (Sikora et 
al. 1994), and it gives a more symmetric Compton flare 
than a — 0. We also obtained fits to the PKS flare with 
a = where the background photons come from the disk 
(Dermer et al. 1992). 

The fit of our model to the PKS 1622-29 data is not 
unique, but the fitting in all cases requires a relatively 
massive black hole (M > 10 9 M©) and a relatively strong 
magnetic field at the base of the jet (B Q > 10 3 G). Also, 
it appears necessary to have an initial ratio of magnetic 
field to rest mass energy-density large compared with 
unity, n ~ 10 - 30. Note that the quasar PKS 1622-297 
is the brightest gamma-ray quasar (Mattox, et al., 1996). 
In other lower gamma luminosity blazars, M and/or B Q 
may be smaller. 

The appearance of one type of flare or another, and also 
the strength of the flares, depends on values of the plasma 
parameters (/ ; °, fi, B Q ) which may be random. Flares 
where the stage of growth is very short may be explained 
by 'spike' type SSC flares (Figure 9a), whereas the slower 
rising gamma light curves may be inverse-Compton flares 
(Figure 9a). Note, that the fast growth and decay of 
luminosity in the beginning of the strong PKS 1622-29 
flare may be connected with 'spike' type SSC flare, which 
appears at smaller initial lepton-ion number ratios (/ ; ° < 
5). 

In the present model, we proposed that the new Poynt- 
ing jet pushes matter continuously. However, if the push 
decreases with time, then the front will lose energy as a 
result of radiation and the luminosity will decrease more 
rapidly than that show here. This effect may also explain 
rapid decay, compared with growth, of the prominent 
flare of the quasar 3C 279 (Kniffen et al. 1993). 

4.2. Radio Variability 

Analysis of the correlation between gamma-ray flares 
and those in radio emissions at short time scales has 
shown that there is no direct correlation in most blazars 
(Reich et al. 1993). Moreover, often the radio flux is 
in the low state during the gamma-ray flare. However, 
observations on longer time-scales show that after strong 
flares in the EGRET band, the radio luminosity starts 
to grow and has a maximum several months later than 
the gamma-ray flare (Reich et al. 1993). This may 
result from either synchrotron self-absorption or free- 
free absorption by plasma surrounding the source core 
(Matveyenko et al. 1992). Thus, the radio emission ap- 
pears at later times when the density in the front and/or 
the external density are much lower. 

In our model the synchrotron self-absorption frequency 
is initially in the infrared or UV, depending on param- 
eters. Thus, initially, there is no radio emission from 
the front. Later, the self-absorption frequency vi(syn) 
slowly decreases (Figure 8) and goes to ~ 3 mm (10 
Hz) after several days. Much later, v\{syn) falls into the 
cm radio band. For some parameters of the model, the 
millimeter emission appears one or two weeks after the 
flare. Initially, only a small fraction of the total luminos- 
ity goes into low frequencies, and the radio flux is very 
small. Because of the slow decrease of v\ (syn) with time, 
vi{syn) ~ 10 12 (i") _1 Hz, the maximum of the radio lu- 
minosity may be a few months later than the gamma-ray 
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Fig. 4. — Total number of leptons Ni and ions JVj in the front and their ratio fu versus time measured by distant observer. 
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Fig. 5. — Number density of leptons n' p ions nj, background photons n' h , synchrotron photons n' syn , synchrotron-self-Compton photons 
n> SSC anc ' mverse -Compton photons n' Com (in units of cm -3 ) in the front frame versus the logarithm of time measured by distant observer. 



flare. Here, we did not take into account possible free- 
free absorption of radio waves by distributed plasma in 
the core of the source exterior to the jet. This absorption 
may also contribute to the delay in the appearance of the 
radio emission. 

The present model of a propagating front in a Poynt- 
ing flux jet is similar in some respects to the model of 
radio variability of Aller, Aller, and Hughes (1985) and 
Hughes, Aller, and Aller (1985) based on weak shocks 
propagating along a hydrodynamic jet. Detailed multi- 
frequency radio observations of the quasar 0528+134 
shows that after a strong gamma-ray flare, the radio flux 
grows and has a maximum few months later (Pohl et al. 
1995). It is important to note that during this event, 
the new VLBI super-luminal component appeared with 
a back-extrapolat-ed time of origin close to the time of 
the gamma-ray flare (Krichbaum et al, 1995). These ob- 
servations support the present model where a strong im- 
pulsive disturbance - the front - propagates along and is 
powered by an otherwise invisible Poynting flux jet. The 



front radiates first in gamma rays, and later becomes 
strong in the radio when it is observed as a new super- 
luminal VLBI jet component. 

5. CONCLUSIONS 

The paper develops a self-consistent dynamical model 
of gamma-ray flares and VLBI outbursts of blazars based 
on lepton acceleration in the propagating front of an oth- 
erwise invisible Poynting flux jet. Inverse-Compton scat- 
tering of the accelerated, relativistic leptons off of back- 
ground photons as well as off of synchrotron photons (the 
SSC process) and synchrotron radiation are taken into 
account. 

It is shown that gamma-ray flares in the frequency 
band of the EGRET instrument (at 50 < E < 10 4 MeV) 
may be of two types: (1.) 'SSC spikes' which have a fast 
rise (~ hr.) and gradual decay of luminosity. This type 
of flare appears in the cases where SSC photons have high 
enough energy to be in the EGRET band. In this case 
inverse-Compton radiation, which usually gives a smaller 
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Fig. 6. — Lorentz factors in the front frame separating different parts of the lepton spectrum: 71,72, and 73 versus time measured by 
distant observer. 7 is the average Lorentz factor. "f r ms is the the root mean square Lorentz factor, "(pairs is the threshold Lorentz factor 
for pair production. For the case shown there is no pair production. 



contribution to the luminosity, is hidden by the strong 
SSC flare. (2.) 'Compton flares' which have a growth 
time of few hours and a duration of days up to weeks 
(depending on the distribution of background photons 
and viewing angle). Comparison with recent observa- 
tions of the gamma-ray flare in PKS 1622-297 (Mattox 
et al. 1996) show that it may be interpreted as a 'Comp- 
ton flare'. To explain the high luminosity and the short 
duration of the PKS flare, the viewing angle must be 
quite small, o b s ~ 0.1 rad. Also, we need a black hole 
mass M~3x 10 9 M Q and a magnetic field at the base 
of the jet B ~ 10 3 G. 

The model gives almost simultaneous flares for photon 
energies from high-energy gamma-ray to lower energy 
gamma, X-ray, UV, visible, and down to the infrared. 
The short wavelength (mm) radio flare for some parame- 
ters may appear shortly (~ days) after the main gamma 
flare, because the self-absorption frequency, which ini- 



tially corresponds to the IR region of spectrum, decreases 
rapidly during a few days (see Figure 8) to short wave- 
length radio band (in some cases, however, mm radiation 
may appear much later). Later, it decreases more slowly, 
so that the centimeter waveband radio radiation in most 
of the cases should appear much later (~ months) than 
the gamma-ray flare. 

In many gamma ray loud objects, both flares and 
steady radiation are observed ( Fichtel and Thompson 
1994; Montigny et al. 1995). The steady radiation may 
come from a superposition of a number of overlapping 
outbursts. In other cases, the EGRET instrument de- 
tects only the highest states of strong flares with the 
weaker steady emission obscured in the noise. 

Extremely high energy photons (~ TeV) have been 
detected from the object Markarian 421 (Punch et al. 
1992). Here, we comment on the possible origin of TeV 
photons in the framework of the present model. TeV 
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Fig. 7. — Total (bolometric) apparent luminosities radiated by the front due to the inverse-Compton, SSC, and synchrotron processes 
(in erg/s) as a function of time measured by distant observer. 



radiation may result from from the SSC process by lep- 
ton scattering off of synchrotron photons if the lepton 
distribution function extends to sufficiently high Lorentz 
factors (~ 10 5 ). However, the magnetic field at the base 
of the jet B D may have to be smaller than the values con- 
sidered here (~ 10 3 G) in order to avoid excessive pair 
production. 
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Fig. 9. — Typical light curves of apparent luminosity in the EGRET band in units of 10 47 erg/s as measured by distant observer. In 
case (a), / ( ° = 10, and the light curve is determined by SSC processes. In case (b), / ; ° = 20, and the light curve is determined by the 
inverse-Compton process. 
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Fig. 10.— Comparison of the EGRET light curve of PKS 1622-297 (the dashed curve) (Mattox ct al. 1996) with the calculated 'Compton 
type' light curve (solid curve). 



